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Abstract 

The reaction mechanism of nucleus- nucleus collisions at projectile 
energies around the Fermi energy is investigated with emphasis on the 
production of fragmentation-like residues. The results of simulations 
are compared to experimental mass distributions of elements with Z 
= 21-29 observed in the reactions 86 Kr+ 124 ' 112 Sn at 25 AMeV. The 
model of incomplete fusion is modified and a component of excitation 
energy of the cold fragment dependent on isospin asymmetry is intro- 
duced. The modifications in the model of incomplete fusion appear 
consistent with both overall model framework and available exper- 
imental data. A prediction is provided for the production of very 
neutron-rich nuclei using a secondary beam of 132 Sn where e.g. the 
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reaction 132 Sn+ 238 U at 28 AMeV appears as a possible alternative to 
the use of fragmentation reactions at higher energies. 

Introduction 

Nucleus-nucleus collisions in the Fermi energy domain exhibit a large variety 
of contributing reaction mechanisms and reaction products ( see e.g. [T] ) 
and offer the principal possibility to produce mid-heavy to heavy neutron- 
rich nuclei in very peripheral collisions. In the reactions of massive heavy 
ions such as 124 Sn+ 124 Sn [2] and 86 Kr+ 124 Sn [3], an enhancement was ob- 
served over the yields expected in cold fragmentation which is at present the 
method of choice to produce neutron-rich nuclei. In this case the neutron- 
rich nuclei are produced in damped symmetric nucleus-nucleus collisions with 
intense nucleon exchange leading to the large width of isotopic distributions. 
Further enhancement of yields of n-rich nuclei was observed in the reaction 
86 Kr+ 64 Ni [1] in the very peripheral collisions, thus pointing to the possible 
importance of neutron and proton density profiles at the projectile and tar- 
get surfaces. In [5], the model of deep-inelastic transfer j6] was supplemented 
with a phenomenological correction introducing the effect of shell structure 
on nuclear periphery. A consistent agreement with experimental data was 
achieved in the reactions of a 25 AMeV 86 Kr beam with three different target 
nuclei, specifically allowing to describe the deviation of the nucleon exchange 
from the path toward isospin equilibration. 

In this article, we present an investigation of the reaction mechanism 
of the nucleus-nucleus collisions with an emphasis on the production of 
fragmentation-like residues at projectile energies around the Fermi energy. 
The fragment at ion- like residues with mass and charge below the projectile 
can be, according to calculation [I], produced at projectile energies around 
the Fermi energy in the process of incomplete fusion. The difference to high- 
energy fragmentation as assumed in standard abrasion-ablation model (0[E] 
) is that here the participant zone fuses with one spectator, thus creating 
a very hot multifragmentation source and the other spectator remains cold. 
The cold projectile-like fragment is not emitted at zero angle but due to clas- 
sical Coulomb trajectory along which the collision evolves it appears at larger 
angles. The existence of such a cold fragment can be deduced e.g. from the 
data of Casini et al. [5] where mass and excitation energy of the projectile- 
like fragment ( PLF ) were kinematically reconstructed and it was established 
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that the products lighter than the projectile are cold and heavier ones are 
hot. Such dependence can not be explained using the concept of nucleon ex- 
change as represented by the model of deep-inelastic transfer [10J. The simu- 
lation using the incomplete-fusion ( ICF ) code |T] reproduces such behavior 
reasonably well, due to low excitation energy of fragmentation-like residues 
with masses below projectile and considerably higher excitation energy of 
ICF-like residues heavier than projectile. Furthermore, a large amount of 
detailed data on hot multifragmentation source in various reactions was re- 
produced, including detailed experimental studies of the mid-velocity source 
[TT] and single fusion-like source [12]. The reaction 124 Sn+ 27 Al was used to 
verify the ICF model using heavy residues with masses A=60-90 measured 
with high precision using a fragment separator at forward angles [T3] and 
both the production cross sections and the N/Z trend were reproduced well 
for these products originating from hot projectile-like ICF source. However, 
little is known on N/Z trend of the cold fragments, which can be produced 
in symmetric reactions of heavy nuclei. In the present paper we examine 
such a trend using the experimentally observed heavy residues in reactions 
86 Kr+ 112 ' 124 Sn at projectile energy 25 AMeV and refine the model description 
in order to improve the prediction of the production of exotic ( neutron-rich 
) nuclei in nucleus-nucleus collisions around Fermi energy. 



Production of heavy residues in reactions 
86 Kr+ ii2,i24 Sn at p ro j ec tile energy 25 AMeV 

According to the incomplete fusion model pQ, in order to study the pro- 
duction of cold fragment at ion- like residues around the Fermi energy it is 
necessary to study symmetric nucleus-nucleus collisions of two massive nu- 
clei ( with target comparable or heavier than projectile ). The products of 
interest will be produced at angles away from zero degrees and the measure- 
ment at these angles is necessary to observe such products. Furthermore, 
according to the calculation, the cold fragment at ion- like residues should be 
increasingly dominant products for the channels with the number of stripped 
protons exceeding seven to eight. Thus the products of primary interest are 
the heavy residues considerably lighter than the initial projectile, which due 
to the removal of a significant number of protons can be also considerably 
neutron-rich. 
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In order to examine the prediction of the ICF code PQ, the results of 
simulations were compared to experimental mass distributions of elements 
with Z = 21 - 29 observed within the separator acceptance in the reaction 
86 Kr+ 124 ' 112 Sn at 25 AMeV [3]. Figs. [HEJshow ^ ne comparison of experimen- 
tal mass distributions ( symbols ) to the results of the simulation ( dashed 
line ) for the reactions of 86 Kr with two tin targets. The simulation uses 
either the model of deep-inelastic transfer [HI [5] for peripheral collisions or 
the model of pre-equilibrium emission and incomplete fusion ( ICF ) for vi- 
olent ( central ) collisions, combined with the de-excitation code SMM [14] . 
The simulated yields were filtered for angular acceptance of the separator 
positioned at 4° with appropriate azimuthal corrections |3j. The descrip- 
tion of experimental data using the de-excitation code SMM, representing 
the combination of statistical multifragmentation model for highly excited 
nuclei with evaporation cascade for lower excitation energies, proved consis- 
tently better than using the model of sequential binary decay. Analogous 
calculation, consisting of the model of deep-inelastic transfer for the early 
stage, combined with the SMM, and thus representing only the peripheral 
collisions, was used for the same reactions in the work [5] and the yields of 
nuclei with Z = 30 - 35 were reproduced well. The contribution from the 
violent collisions is negligible for nuclei with Z=30-35 and thus description of 
violent collisions used in the present work does not affect the results for these 
reaction products. It is however essential for description of the production of 
nuclei around nickel and below where the ICF contribution becomes increas- 
ingly dominant. One can see that for the reaction 86 Kr+ 124 Sn the calculation 
overpredicts the experimental yields of the most neutron-rich products below 
nickel ( with more than 8 stripped protons ). The situation is similar also 
for the reaction 86 Kr+ 112 Sn even if the most neutron-rich products below 
nickel are less populated than in reaction of 86 Kr+ 124 Sn. Thus the yields 
of neutron-rich nuclei appear to be overpredicted in the region where the 
cold fragment at ion- like residues dominate. A modification of the model of 
incomplete fusion concerning the cold fragmentation-like residues should be 
considered, specifically concerning the description of the N/Z-trend which 
was not verified in the past. 

The model of incomplete fusion p], used in the calculation, considers a 
spectator-participant scenario evolving along the classical Coulomb trajec- 
tory, followed by fusion of the participant zone with one spectator, typically 
the heavier one due to larger contact area and thus larger attractive force. 
The charge of spectator zones is determined using the combinatorial probabil- 
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ity, which is a standard approach in fragmentation codes [SJ. The excitation 
energy of the cold fragment is determined considering the two-body colli- 
sions of participant and spectator nucleons along the separation plane pQ. 
The concept of combinatorial probability explores the available statistical 
phase space and allows a rather wide range of isospin asymmetries. From a 
dynamical point of view, however, the spectator-participant scenario implies 
an instantaneous separation and thus preservation of isospin asymmetry of 
the initial nucleus in the ground state, with homogeneous density in the in- 
terior of the nucleus. Thus, the two concepts seem to be in contradiction, 
which can be resolved when assuming that the change of isospin asymmetry is 
dynamically consistent with transfer of certain amount of nucleons across the 
separation plane. Simultaneously, the relative velocity between the partici- 
pant and spectator zone increases from zero to maximum value corresponding 
to final incomplete-fusion scenario and thus the transferred nucleons should 
carry this relative velocity which will be transferred into excitation energy of 
the acceptor. Thus a component of excitation energy dependent on isospin 
asymmetry can be deduced. It can be assumed that, due to absence of 
Coulomb barrier, the nucleons transferred will be predominantly neutrons. 
The formula for such isospin dependent component of excitation energy of a 
cold spectator ( acceptor ) can be written as 

V ICF rp _y 

E S (A S , Z 3 ) =x{A s - A (Zs)) (-^) 2 ^— 1° (1) 

"proj **-P 

where Ep, Ap are the projectile kinetic energy and mass, Vc is the 
Coulomb barrier, v pro j, v^ F are the projectile velocity and the final relative 
velocity between hot and cold fragment in the incomplete-fusion scenario, 
evaluated at the Coulomb barrier, Ag, Zg are mass and charge of the spec- 
tator ( cold fragment ), A Q (Zg) is the spectator mass corresponding to N/Z 
of initial nucleus and a; is a random number between zero and one, gener- 
ated for each collision. The random number is introduced due to uncertainty 
concerning the exact moment of transfer and represents an zero-th order es- 
timate allowing to reproduce the mean value of extra excitation energy due 
to transfer of neutrons. The excitation energy is evaluated only in the case 
when the cold fragment is acceptor ( and thus neutron- rich ) . For the loss of 
neutrons no such component is evaluated, thus assuming that neutrons close 
to the Fermi level are lost and no enhancement of intrinsic excitation occurs. 

The results of modified calculation employing the formula ([I]) are shown 
in Figs, dj [2] as solid lines. One can see that the overall agreement with the 
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experimental data is improved in the modified calculation compared to the 
standard one. The calculation with modification, which is fully consistent 
with overall description of the nucleus-nucleus collisions and introduces no 
free parameters, describes the trend of the experimental yields of the most 
neutron-rich products below nickel for the reaction 86 Kr+ 124 Sn, which the 
initial calculation overpredicts. The situation improves also for the reaction 
86 Kr+ 112 Sn. For the most neutron-rich products below chromium nuclei oc- 
casionally appear less populated than it is predicted by the model. This 
however can be caused by missing yield in the experiment due to high back- 
ground from the initial beam, because of increasing overlap of the charge 
states ( in terms of magnetic rigidity ) with charge states of the scattered 
beam. The effect is more pronounced for the reaction 86 Kr+ 112 Sn due to 
lower experimental yields. 

It is worthwhile to mention that both the original and modified calcula- 
tions correctly reproduce the experimental cross sections of heavy residues 
with Z = 25 - 30 for the reaction of 86 Kr+ 64 Ni at 25 A MeV @], measured at 
angles below 3 °, as documented in Fig. [3] ( dashed and solid line represent 
original and modified calculation ). The discrepancies at the proton-rich side 
are caused by restricted Bp-coverage for such products in the experiment 
focused primarily on neutron-rich products. The situation in Fig. [3] and its 
comparison with reactions 86 Kr+ 112 ' 124 Sn shows that indeed measurement at 
larger angles and use of heavier target is necessary in order to observe the 
experimental cross section of cold fragment at ion- like residues. 

In general, the modifications in the model of incomplete fusion appear 
consistent with both overall model framework and experimental data and 
thus one can expect improved predictive power which can be used to predict 
production of exotic mid-heavy to heavy neutron-rich nuclei in the reactions 
around the Fermi energy, and possibly identify under which conditions such 
approach can be more effective than other methods. From the point of view 
of reaction dynamics, the modified model of incomplete fusion is consistent 
with the formation of a neutron-rich region between cold and hot fragment 
( or participant zone as its precursor ). Similar effect was reported in the 
literature [TJ] as a possible consequence of the evolution of nuclear mean field. 
The number of transferred neutrons can then be determined by a mechanism 
similar to the random neck rupture, as established in nuclear fission [16J, 
which can justify the applicability of a combinatorial ( and thus essentially 
statistical ) probability in the description of dynamical reaction mechanism 
such as the incomplete fusion. 
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Production of neutron-rich nuclei around N=82 



A great deal of attention was paid in recent years to production of the sec- 
ondary beams of exotic nuclei. One of the most promising ways to produce 
extremely neutron-rich nuclei around the neutron shell N=82 is fragmenta- 
tion of a secondary beam of 132 Sn. Nevertheless, based on the results of the 
previous section, one can in principle consider also the reaction in the Fermi- 
energy domain at energies below 50 AMeV. The comparison of production 
cross sections for the reaction 132 Sn+ 238 U at 28 AMeV with fragmentation 
cross section of 132 Sn beam with Be target is provided in Fig. HI For the reac- 
tion 132 Sn+ 238 U the modified DIT code [5] was used for peripheral collisions 
together with original model of incomplete fusion pQ ( dashed lines ) and 
its modification presented in this article ( solid lines ) for central collisions, 
while for the fragmentation of 132 Sn beam the codes COFRA [17] ( dotted 
lines ) and EPAX [18] ( dash-dotted lines ) were used. The production cross 
sections calculated using both the original and modified model of incomplete 
fusion for Z=46 are comparable with results of EPAX and COFRA, while 
for elements with lower atomic numbers the reaction 132 Sn+ 238 U leads, ac- 
cording to both the original and modified model of incomplete fusion, to still 
more favorable cross sections exceeding both COFRA and even EPAX cross 
sections. The high cross sections for proton-stripping channels with 238 U tar- 
get were observed recently in reaction 58 Ni+ 238 U by Corradi et al. [TO] at 
energies about 1 AMeV above Coulomb barrier. The proton-stripping cross 
sections were reproduced using the calculation analogous to the one used here 
after minor readjustment ( increase by 0.5 fm ) of diffuseness reflecting mean- 
field effects at such low energy [20] • When looking at the calculated cross 
sections for 132 Sn+ 238 U in Fig. H] one can see that the trend of the shapes 
of mass-distributions calculated using original model of incomplete fusion is 
quite similar to EPAX ( since the assumptions are essentially similar ) while 
trend of the elemental yields differs due to assumed classical Coulomb trajec- 
tories vs straight movement assumed in fragmentation model. The validity of 
EPAX code concerning the trend of n-rich products is a matter of discussion 
and based on the results of the present work one can expect that experimental 
production cross sections of very n-rich nuclei in fragmentation reactions will 
be significantly lower than the prediction of EPAX parametrization. The 
magnitude of such effect needs to be verified ( a dedicated experiment on 
fragmentation of 132 Sn secondary beam is planned at GSI [21] ), one can 
nevertheless expect that the code with more realistic treatment of specta- 
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tor excitation energy, such as COFRA, will provide better predictions. It is 
interesting to note that the difference of predictions of fragmentation cross 
sections by EPAX and COFRA is similar to difference of predictions of pro- 
duction cross sections in reaction 132 Sn+ 238 U by the original and modified 
model of incomplete fusion. Despite different processes, the amount of exci- 
tation energy of the spectator appears to be comparable in both cases and 
the success of modified model of incomplete fusion in the present work may 
imply that COFRA will be successful in the fragmentation reactions. 

The in-target yields calculated using the production cross sections from 
Fig. H]are shown in Fig. [5j For the fragmentation of 132 Sn secondary beam 
an energy 100 AMeV was used, which is the maximum foreseen for post- 
accelerator envisioned for future European RNB facility Eurisol [22]. The 
achievable in-target reaction rate was determined using code AMADEUS 
[23] . For the reaction 132 Sn+ 238 U at 28 AMeV a target thickness 40 mg/cm 2 
was assumed. For the intensity of 132 Sn secondary beam a value of 10 12 s _1 
was adopted from Eurisol RTD Report [22]- Due to larger target thickness, 
the in-target yield for fragmentation option calculated using both COFRA 
and EPAX dominate for elements Z=44 and above, for lighter nuclei never- 
theless the larger production cross sections in the Fermi-energy domain lead 
also to larger in-target yields despite relatively thin target and for Z=40 the 
in-target yield calculated using the modified model of incomplete fusion ex- 
ceeds the COFRA value ( and the EPAX value is exceeded by original model 
of incomplete fusion ). To answer on viability of Fermi-energy option, a cru- 
cial question is which calculation of fragmentation cross sections is realistic. 
Recent measurement carried out at MSU [21] has shown that experimental 
fragmentation cross section of neutron-rich Ni isotopes are overpredicted by 
EPAX by up to two orders of magnitude. If the situation is analogous in 
the fragmentation of 132 Sn the Fermi-energy option can become interesting. 
However, the angular distribution of reaction products at 28 AMeV would 
require a large-acceptance separator with angular coverage up to 10 degrees 
and a highly efficient gas-cell in order to form a secondary beam. 

Conclusions 

The reaction mechanism of the nucleus-nucleus collisions at projectile ener- 
gies around the Fermi energy was investigated with emphasis on the produc- 
tion of fragment at ion- like residues. The results of simulations were compared 
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to experimental mass distributions of elements with Z = 21 - 29 observed 
around 4° in the reaction 86 Kr+ 124 ' 112 Sn at 25 AMeV. The model of incom- 
plete fusion was modified and a component of excitation energy of the cold 
fragment dependent on isospin asymmetry was introduced. The modifica- 
tions in the model of incomplete fusion appear consistent with both overall 
model framework and available experimental data. Prediction is provided 
for the production of very neutron-rich nuclei using a secondary beam of 
132 Sn where e.g. the reaction 132 Sn+ 238 U at 28 AMeV appears as a possible 
alternative to the use of fragmentation reactions at higher energies. 
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Figure 1: Comparison of the simulations to experimental mass distributions (symbols ) 
of elements with Z = 21 - 29 observed around 4° in the reaction 86 Kr+ 124 Sn at 25 AMeV 
[3j . Dashed line - results of the standard simulation [IJ [5] combined with the de-excitation 
code SMM [14] , Solid line - results of simulation using modified model of incomplete fusion 
(eq. ffli i. 
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Figure 2: Comparison of the simulations to experimental mass distributions (symbols ) 
of elements with Z = 21 - 29 observed around 4° in the reaction 86 Kr+ 112 Sn at 25 AMeV 
[3]. Solid, dashed lines - as in Fig. [1] 



12 







Z=30 

V. 








Z=29 




Z = 28 


g 5 




/ \" 
/ " \* 

If \ 


(mb) 

0.1 1 






/" \ 


«a s 


/ 

/ ■ ■! 


b 






b „ 






it \ 


b H 






















50 


60 70 SO 

A 




50 




60 70 80 

A 




50 60 70 80 

A 






Z=27 








Z = 26 




Z = 25 


i 2 




■ \. 


(mb) 

0.1 1 








(mb) 

0.1 1 


' / 


b _ 


- i 


* \" 
1 . 

» 

V 


b _ 






T 

V 


b H 


ft \ 




50 


60 70 ao 




50 




60 70 80 




50 60 70 80 



AAA 

Figure 3: Comparison of the simulated ( lines ) and experimental ( symbols ) cross 
sections of heavy residues with Z = 25 - 30 for the reaction of 86 Kr+ 64 Ni at 25 A MeV 
[4], measured at angles below 3 °. Dashed and solid line represent standard and modified 
calculation. 
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Figure 4: Comparison of production cross sections for reaction 132 Sn+ 238 U at 28 AMeV 
using standard ( dashed lines ) and modified simulation ( solid lines ) with fragmentation 
cross sections of 132 Sn beam with Be target using COFRA [T7] ( dotted lines ) and EPAX 
[18] ( dash-dotted lines ). 
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